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(N ■ Abstract: 

Oh 1 We study finite volume effects in heavy quark systems in the framework of heavy meson chiral perturbation theory for 
full, quenched, and partially quenched QCD. A novel feature of this investigation is the role played by the scales A* 
and 5 S , where A* is the mass difference between the heavy- light vector and pseudoscalar mesons of the same quark 
' content, and 5 S is the mass difference due to light flavour SU(3) breaking. The primary conclusion of this work is 
that finite volume effects arising from the propagation of Goldstone particles in the effective theory can be altered 
fSJ ' by the presence of these scales. Since A* varies significantly with the heavy quark mass, these volume effects can 
^ \ be amplified in both heavy and light quark mass extrapolations (interpolations). As an explicit example, we present 
results for B parameters of neutral B meson mixing matrix elements and heavy-light decay constants to one-loop order 
in finite volume heavy meson chiral perturbation theory for full, quenched, and Nf = 2 + 1 partially quenched QCD. 
. Our calculation shows that for high-precision determinations of the phenomcnologically interesting SU (3) breaking 
' ratios, finite volume effects are significant in quenched and not negligible in partially quenched QCD, although they 
■^j- . are generally small in full QCD. 

o 
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I. INTRODUCTION 

Numerical calculations of hadronic properties using lattice QCD have provided significant inputs to particle physics 
phenomenology. In particular, the joint effort between experiment and theory to investigate the unitarity triangle in 
the Cabibbo-Kobayashi-Maskawa (CKM) matrix from B meson decays and mixing has made impressive progress 
in which lattice QCD has played an important role. Nevertheless, current lattice calculations are still subject to 
various systematic errors. In this paper, we address finite volume effects which arise in lattice calculations for heavy- 
light meson systems from the light degrees of freedom. Our framework is heavy meson chiral perturbation theory 
(HMxPT) with first order l/Mp and chiral corrections. We assume the mass hierarchy 

A/ GP < A x < M P , (I) 

where Mgp is the mass of any Goldstone particle, Mp is the mass of the heavy-light meson, and A x is the chiral 
symmetry breaking scale. Under this assumption, we discard corrections of the size 

(2) 

M P V ; 

Concerning the finite volume, we work with the condition that 

M GP L » I, (3) 

where L is the spatial extent of the cubic box. Therefore, given that UL/V2 (f n « 132 MeV) will be close to one m 
lattice simulations in the near future, one can still neglect the chiral symmetry restoration effects resulting from the 
Goldstone zero momentum modes 0, E| when Eq. © is satisfied. 

The main task of this work is to study the volume effects due to the presence of the scales 

A, = M P * - M P , (4) 

and 

S s = M Ps - M P , (5) 

where P* and P are the heavy- light vector and pseudoscalar mesons containing auorrf anti-quark 1 , and P s is the 
heavy-light pseudoscalar meson with an s anti-quark. The scale A» appears due to the breaking of heavy quark spin 
symmetry that is of 0(\/Mp) and S s comes from light flavour SU(3) breaking in the heavy-light meson masses. Under 
the assumption of Eq. JQ), A* is independent of the light quark mass, and S s does not contain any l/Mp corrections, 
at the order we are working. 

In the real world, both A* and S s are not very different from the pion mass. In fact Q, 

M Bs - M B = 91 MeV, (6) 
M Ds - M D = 104 MeV, (7) 
M B > - M B = 46 MeV, (8) 



and 



M d* - M D = 142 MeV. (9) 



In current lattice simulations, these mass splittings vary between and ^150 MeV. Therefore it is important to 
include them in the investigation of finite volume effects. Equation @ implies that the Compton wavelength of the 
Goldstone particle is small compared to the size of the box. Therefore finite volume effects mainly result from the 
propagation of the Goldstone particles to the boundary. However, as shown in Section IIIII A* and 8 S can, in a 



1 We work in the isospin limit in this paper. 
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non-trivial way, alter these effects. In particular, since A* varies with the heavy quark mass, finite volume effects can 
be significantly amplified in heavy quark mass extrapolations. 

This paper is organised as follows. In Section |n] we summarise the ingredients of HMxPT relevant to this work. 
Section lTTTl is devoted to the discussion of HMxPT in finite volume, emphasising the role played by S s and A* . We then 
present an explicit calculation of neutral B meson mixing and heavy- light decay constants in Section IIVI and discuss 
the phenomenological impact that finite volume effects can have. We conclude in Section E] Some mathematical 
formulae and results are summarised in the appendices. 

As this work progressed, we were informed that similar ideas and techniques were also being applied in heavy baryon 
chiral perturbation theory [f| 2 . Although the underlying physics is somewhat different, many technical aspects are 
quite similar to those presented here. 



II. HEAVY MESON CHIRAL PERTURBATION THEORY 

The chiral Lagrangian for the Goldstone particles is 

£gp = ^(s)tr[(^Et)(^E)+Et x + x t E ] 

+A V , [a(d^ )(d^o) - Mq^o] . (10) 

where A n > = 1 for quenched QCD (QQCD) and partially quenched QCD (PQQCD), and A v > = for full QCD. 
E = exp(2i<f>//) is the non-linear Goldstone particle field, with $ being the matrix containing the standard Goldstone 
fields. 3 We use / = 132 MeV. In this work, we follow the supersymmetric formulation of (partially) quenched chiral 
perturbation theory [(P)QxPT] 0,0. Therefore E transforms linearly under SU(3) L ® SU(3) R , SU(3\3) h ®SU{3\3) R 
and SU(6\3) L <g> SU(6\3) R in full QCD, QQCD and PQQCD respectively. The symbol "(s)tr" in the above equation 
means "trace" in full QCD and "supertrace" in (P)QQCD. The variable x is defined as 

n „ , , -2(Q\uu + ddlO) , , 
X ee 2B M g = V 1 —M q , (11) 

where the quark mass matrix M. q is 



in full QCD, 



^(QCD) = diag ( TO; 772, m s ), (12) 
A'j (QQCD) = diag(m, m, m s ,m, m, m s ), (13) 

valence ghost 

in QQCD, and 

M. ( p QQ CD ) = diag(m, m, m s , rh, rh, fh s , m, m, m s ), (14) 

valence sea ghost 

in PQQCD. We keep the strange quark mass different from that of the up and down quarks in the valence, sea and 
ghost sectors. Notice that the flavour singlet state $o = str($)/V6 is rendered heavy by the U{\)a anomaly in 
PQQCD HG3 

and can be integrated out; it has to be kept as a dynamical degree of freedom in QQCD. 
The inclusion of the heavy- light mesons in chiral perturbation theory (HMxPT) was first proposed in Refs. 0, 
Il2l Il3j. with the generalisation to quenched and partially quenched theories given in Refs. 0,^3- The 1/Mp and 
chiral corrections were studied by Boyd and Grinstein |ig in full QCD and by Booth [l7| in QQCD. The spinor field 
appearing in this effective theory is 

H (Q) = 111 



(P&>-f ~ ^ (Q) 7 5 ) , (15) 



2 We thank Silas Beane for drawing our attention to his work. 

3 In this paper, we only address situations where there are no multi-particle thresholds involved in loops. This is the case for the explicit 
calculation presented in Section ll VI Therefore, in spite of the sickness pointed out in Ref. 6i|, we can still use the Minkowski formalism 
even for the case of (P)QQCD. This makes the physics discussion in Section IlIII simpler compared to the Euclidean formalism. The 
effects from multi-particle thresholds in finite volume HMxPT will be discussed in a future publication. 
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where pP and Paffi annihilate pseudoscalar and vector mesons containing a heavy quark Q and a light anti-quark 
of flavour a. Under a heavy quark spin SU (2) transformation S, 

Hjf»^SH<?». (16) 

Under the vector light-flavour transformation U [i.e., U € SU(3) for full QCD, U G SU(3\3) for QQCD and U E 
SU (6\3) for PQQCD], 

H (Q)^ H (Q)UI (17) 

Also, the conjugate field, which creates heavy-light mesons containing a heavy quark Q and a light anti-quark of 
flavour a, is defined as 

Hp =7 o H (Q) 10 . (18) 
Furthermore, the Goldstone particles appear in the HM^PT Lagrangian via the field 

t = e**", (19) 

which transforms as 

£ — ► l^t/ 1 = (20) 

where C^l(r) is an element of the left-handed (right-handed) SU(3), SU(3\3) and SU(6\3) groups for QCD, QQCD, 
and PQQCD respectively. The HM^PT Lagrangian, to lowest order in the chiral and 1/Mq expansion, for mesons 
containing a heavy quark Q and a light anti-quark of flavour a is then 

£hm x pt = -ztr D (i^)^i^ 

+ i ff tr D (SP W [^d^-td^^H^+B^^tm (h^hPws) str - $d^] , 

(21) 

where B n ' = for full QCD, and B n > — 1 for (P)QQCD 4 . We do not distinguish the coupling g in these theories. It is 
implicitly assumed that "(s)tr Q " is taken appropriately in flavour space. trD means taking the trace in Dirac space. 
The HM%PT Lagrangian for mesons containing a heavy anti-quark Q and a light quark of flavour a is obtained by 
applying the charge conjugation operation to the above Lagrangian [lc|. At this order, the propagators for Pjp^ and 

r>*(Q) 

P a mesons are 



2(vk + ie) ' 2(v-k + ie) ' 
respectively. 

The effects of chiral and heavy quark symmetry breaking have been systematically studied in full 0] and 
quenched 17] HMxPT. Amongst them, the only relevant feature necessary for the purpose of this work, i.e., the 
investigation of finite volume effects, are the shifts to the masses of the heavy-light mesons. These shifts are from the 
heavy quark spin breaking term 



A 2 



M P 

and the chiral symmetry breaking terms 

\ X B Q tr D (BP [tMgt + tfMg?]^ 
+X[B tr D (B^hP) [W 9 t + ?M q t% 



tr D (B^a^H^a^) , (23) 



r(Q) 



(24) 



4 However, since we integrate out the rf in PQQCD, the coupling 7 does not appear in the results presented in this paper. 
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We choose to work with the effective theory in which the heavy-light pseudoscalar mesons that contain a heavy quark 
and auord valence anti-quark are massless. Notice that the term proportional to X[ in Eq. I|24ll causes a universal 
shift to all the heavy-light meson masses. This means that the masses appearing in the propagators of heavy vector 
mesons and any meson containing an s anti-quark (valence or ghost) arc shifted as follows: 

-i{g^ ~ v M v u ) i 



2(v ■ k - A* + if) ' 2(vk-5 s +ie)' 
and 

2(v ■ k - A* - S s + ie)' 



(26) 



for P* , P s , and P* (heavy vector meson containing an s valence or ghost anti-quark), respectively. The mass shifts 
can be written in terms of the couplings in Eqs. i|23[l and (|24|l : 

A, - -8*, (27) 

and 

S s = 2\ 1 B (m a -m). (28) 

In PQQCD, there are two additional mass shifts because the sea quarks have different masses from the valence and 
ghost quarks: 

5 s = Mp a -Mp = 2\iB (m s -m), (29) 

and 

S sea = Mp-M P = 2X 1 B (m-m). (30) 

where P (P s ) is the heavy- light pseudoscalar meson with a d (s) sea anti-quark. The propagators of the heavy mesons 
containing sea anti-quarks are: 

(31) 



2(vk — 4ea + *e) ' 

—i{9nv - VpVv) 



(32) 



2(v -k — S sca - S s + if) 
and 



2(v ■ k - A* - S sca + if) ' 

(33) 



- (34) 
2(v ■ k - A, - 5 sea - S s + if) 

for P, P* (vector meson with a d sea anti-quark), P 8 , and P* (vector meson with an s sea anti-quark), respectively. 



III. FINITE VOLUME EFFECTS 



In this section, we discuss generic features of finite volume effects in HMxPT. For clarity, we use the symbol A for 
one of (A*, S s , 5 a , S sca ) or any sum amongst them. 

In the limit where the heavy quark mass goes to infinity and the light quark masses are equal, all the heavy mesons 
in HM^PT become on-shell static sources, and there is a velocity superselection rule when the momentum transfer 
involved in the scattering of the heavy meson system is fixed |l9| . For illustration, consider the vertex with coupling 



6 



g in £hmxPT introduced in Eq. i|21[l . The heavy- light meson P can scatter into by emitting a Goldstone particle 
with mass Mqp through this vertex. The momenta of the mesons P and P* s ■. are 

M pVfl , (35) 

and 

M P * e) w M + fc p = Mpiip + (36) 

where the velocity = (1,0,0,0) in the rest frame of the heavy mesons, and k^ is the soft momentum carried by 
the Goldstone particle. The infinitely heavy P and Pf-, mesons do not propagate in space. Therefore, when such a 
system is in a cubic spatial box, finite volume effects result entirely from the propagation of the Goldstone particle 
to the boundary with momentum k ~ Mqp. In this case, the volume effects behave like exp(— MqpL) multiplied by 
a polynomial in l/L. 

The breaking of heavy quark spin and SU(3) light flavour symmetries in HMxPT can induce a mass difference 

M P * = M P + A, (37) 

which complicates the above picture. In this scenario, the P*^ is still regarded as a static source, but it is off-shell 
with the virtuality A. The period during which the Goldstone particle can propagate to the boundary is limited by 
the time uncertainty conjugate to this virtuality, i.e., 

(38) 

This means that finite volume effects, which arise from the propagation of the Goldstone particles in such a system, 
will decrease as A increases. Equation also indicates that the suppression of the volume effects by a non-zero A 
is controlled by the parameter 

(39) 

To see explicitly how this phenomenon appears in a calculation, we consider a typical sum in one-loop HMxPT, 
with a Goldstone propagator and a heavy-light vector meson propagator in the loop, in a cubic box with periodic 
boundary condition: 

J(M Gr ,A) - (4„) 

k 

where the spatial momentum k is quantised in finite volume as 

(4i, 

with i being a three dimensional integer vector. Using the Poisson summation formula, it is straightforward to show 
that 

J{M GP , A) = J(Af GP , A) + J fv (Mgp, A), (42) 

where 

J (Mqp, A) = -iJ (fc2 _ M ^ + ie)(v ■ k - A + ie) ' (43) 



is the infinite volume limit of J^(Mqp, A), and 



n#0 
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with 



= \l \k\ 



(45) 



is the finite volume correction to J(Mqp, A). In the asymptotic limit where MqpL 1 it can be shown that (with 
n = \n\) 



where 
A 



e [z > [1 - Erf (z) ] 



1 



uMqpL 



J FV (M GP ,A) 



1 



nMspL 
9z 5z 3 



V?r V64 32 



1 (z z 3 



i 



\ 8nnL 



(46) 



^ 2 )[1-Erf(z)] 





1 




+ ( 


1 




uMqpL 





,(47) 



with 



/ A 



GP 



nM GP L 



(48) 



The quantity A is the alteration of finite volume effects due to the presence of a non-zero A. It multiplies the 
factor exp(— uMqpL), which results from the propagation of the Goldstone particle to the boundary. It is possible 
to analytically compute the higher order corrections of A in powers of 1/ ' (uMqpL). This way, one can achieve any 
desired numerical precision. Here it is clear that this alteration of volume effects is controlled by the ratio in Eq. 139fl . 

Next, we consider different limits of A at fixed M G p and L. When A = 0, clearly A = 1. If A is very small 
compared to Mgp, such that z -C 1, A is dominated by the (1/MqpL) term, i.e., 



A m e {z ^ [1 - Erf (2)] 



(49) 



Since Erf(z) grows much faster than exp(z 2 ) in this regime, A will decrease as A increases. When A is of O(Mqp) 
or larger, z » 1, and one can perform an asymptotic expansion of the error function. It can be shown that in this 
situation, 



z 



(50) 



That is, A also decreases as A increases. We have also numerically checked that this is true when z « 1. This means 
that the asymptotic formula in Eq. 1)46(1 reproduces the physical picture outlined in the beginning of this section for 
any A. To demonstrate how fast the asymptotic form in Eq. 147(1 converges to Eq. 144(1 . we define 



dJ F v(M G P,A) 



JfTWgp,A) 



Tasymp / 



GP 



,A) 



J^ m (M GP ,A) 



(51) 



where J^ m 



is the function Jpv evaluated numerically [Eq. l(4~4T) ]. and Jpy mp is the asymptotic form in Eq. (|4*7jl . In 



Fig.^ we plot dJpy as a function of Mqp with three choices of A. It is clear from this plot that Jfv is approximated 
well (to < 3%) by the asymptotic form when MqpL > 2.5. We use the asymptotic forms for integrals of this type 
throughout this work. Also, in this paper we only include the terms with \n\ = 1, y/2, y/3, ^fi and -s/5 in the Poisson 
summation formula. We have confirmed that truncating the sum at \n\ = \fS is a very good approximation (to ~ 3%) 
when MqpL > 2.5. The function Jfy(Mqp, A) is plotted against Mqp in Fig. [3 with L = 2.5 fm and three choices 
of A. It is clear from this plot that A can significantly alter the finite volume effects in J{Mqp : A). 
Another typical sum that appears in one-loop HMxPT in finite volume is 



/C(M GP ,A) 



E 



dko 
~2T 



k - A + ie 



(52) 



It is straightforward to show that 



JC(M G p, A) = K(M G p, A) + Kpv (Mqp , A), 



(53) 
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dJ FV (M GP ,A) 




FIG. 1: d Jfv (Mgp, A), defined in Eq. (1511 . plotted against Mgp, with three different choices of A. This function indicates 
the deviation (in percent) of the asymptotic form of Jfv from the definition in Eq. 1441 . The size of the volume in this plot 
is L =2.5 fm. The Goldstone mass Mgp = 0.197 GeV corresponds to MgpL = 2.5, and Mgp = 0.32 GeV corresponds to 
M GP L = 4. The curve for A = 0.1 GeV is hidden by that for A = 0.2 GeV. 



J FV (M GP ,A) 



0.0025 


\ 


0.002 


\ 
\ 

A 


0.0015 


0.001 




0.0005 


\> * 



A = 0.00 GeV 

A = 0.10 GeV 

A = 0.20 GeV 




M GP (GeV) 



FIG. 2: Jfv (Mgp, A) plotted as a function of Mgp, with three different choices of A corresponding to the three curves. The size 
of the volume in this plot is L =2.5 fm. The Goldstone mass Mgp = 0.197 GeV corresponds to MgpL = 2.5, and Mgp = 0.32 
GeV corresponds to MgpL = 4. 



where 



K(MoP ,A) = ffig^> (54) 

is the infinite volume limit of /C(Mqp, A) and 

<9 Jfv (Mgp, A) 

K FV (M GP ,A) = — (55) 

is the finite volume correction to K(Mqp,A). The function K F y(Mqp, A) is plotted against Mqp in Fig. [3] with 
L = 2.5 fm and three choices of A. As expected, \K fv (Mq P , A)| also decreases when A increases at fixed M G p and 
L. 



IV. NEUTRAL B MIXING AND HEAVY-LIGHT DECAY CONSTANTS 



The study of neutral B meson mixing allows the extraction of the magnitude of the CKM matrix element Vtd, 
and hence the determination of one of the sides of the unitarity triangle. The frequency of the Bd—Bd oscillations, 
which is given by the mass difference, Am^, in this mixing system has been well measured by various experimental 
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FIG. 3: Kfv(Mgp, A) plotted as a function of Mgp, with three different choices of A corresponding to the three curves. The 
size of the volume in this plot is L = 2.5 fm. The Goldstone mass Mgp = 0.197 GeV corresponds to MgpL — 2.5, and 
Mgp = 0.32 GeV corresponds to MgpL = 4 in this plot. 



collaborations It is also calculable in the standard model via an operator product expansion in which the top 
quark and W boson are integrated out. Resumming the next-to-leading order (NLO) short-distance QCD effects, one 
obtains 

GF ;M^\Vt d V t l\ 2 7 lB So(xt)CB(p) 



Am d 



8tt 2 



\(B d \Of B = 



\^)\B d )\ 



2M F 



(56) 



where [i is the renormalisation scale, Xt = m 2 /Mfr, and S (x t ) « 0.784x°- 76 (to better than 1%) is the relevant 
Inami-Lim function j2jj. The coefficients tjb = 0.55 and Cb(jj) are from short-distance QCD effects 0,H3|- The 
matrix element of the four-quark operator 

0AB=2 = [h , {1 _ 7g)d][ 5 7M (l _ 75 ) d ] (57) 

between B d and B d states contains all the long-distance QCD effects in Eq. (|56|) . and has to be calculated non- 
perturbatively. Since |Vtb| = 1 to good accuracy and Am d has been well measured, a high-precision calculation of 
{B d \0^ B=2 {y)\B d ) enables a clean determination of \Vtd\- 

The frequency of the rapid B s —B s oscillations can be precisely measured at the Tevatron and LHC Qj. Therefore 
an alternative approach is to consider the ratio 

1 ' " ' (B s \Of B = 2 \B s 



Am s 
Am d 



Vtf 



v. 



ill 



Me 
Mr 



(B d \0* 



AB=2 



B d 



(58) 



in which many theoretical uncertainties cancel. Here Am s is the mass difference in the B s —B s system and Of B=2 — 
[67^(1 — 75 )s] [&7 M (1 — 75)5]. The unitarity of the CKM matrix implies |Vt s | ~ \V c b\ to a few percent, and \V c b\ can be 
precisely extracted by analysing semileptonic B decays Therefore a clean measurement of Am s /Am d will yield 
an accurate determination of \V fd \. 

The matrix elements in Eq. I|58|l are conventionally parameterised as 



{B q \Of s = 2 



B q ) = ^bJbBbM' 



(59) 



where the parameter Bs q measures the deviation from the vacuum-saturation approximation of the matrix element, 
and q — d or s. The decay constant fs q is defined by 

(0\h fll5 q\B q (p))=ip fl f Bq . (60) 

Lattice QCD provides a reliable tool for calculating these non-perturbative QCD quantities from first principles 5 . 



5 Some selec ted reviews in the long history of lattice calculations for the B mixing system can be found in Refs. |2.'l I24L l2fl 12(1 l27t I2SL 
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Since Am s / Arrid will be measured to very good accuracy, it is important to have clean theoretical calculations for 
[the SU(3) breaking ratios of] the matrix elements, decay constants and B parameters involved. Current lattice 
calculations have to be combined with effective theories in order to obtain these matrix elements at the physical quark 
masses. This procedure can introduce significant systematic errors and dominate the uncertainties in the SU (3) 
breaking ratio j3j> |33| 



e= /s 1 V^ ) (61) 



which is the key theoretical input for future high-precision determination of \Vtd\ via the study of neutral B mixing 6 . 
Howe ver, the use of effective theory also offers a framework for studying finite volume effects in lattice calculations [a, 
IE IM HE HE H3, US E3, E3, E2, E3, 13, EE E^. We will demonstrate in this section that finite volume effects 
might turn out to exceed the current quoted systematic errors for quantities such as £. 

A. The one-loop calculation in finite volume 

In this subsection, we discuss one-loop calculations for the B parameters and heavy-light decay constants mentioned 
above in finite volume HM%PT including the appropriate mass shifts to the first non-trivial order of the chiral and 
1/Mp expansion. The inclusion of other first-order corrections in these quantities is straightforward. It simply 
introduces additional low-energy constants (LECs) which account for short-distance physics and do not give rise to 
finite volume effects at this order, so we will not discuss this issue here. We have performed the calculation for full 
QCD, QQCD and PQQCD with the mass shifts given between Eqs. (J25} and 

For the purpose of this work, the axial current bj^j^qa is 



A, 



(| ) tr D [wX Q) £] , (62) 



and the four-quark operator O ap «- 2 (when P a = B dmS , Q AP *- 2 becomes Off ~ 2 ) is 



Qaa 



4/3[(ep; (Q)t ) a (ep* (( ^ 



(63) 



in HM^PT [l8| , where k and /3 are the low-energy constants which have to be determined from experiments or lattice 
calculations. Notice that the index a in Eq. I|63|) is not summed over. Again, the inclusion of the chiral and 1/Mp 
corrections in these operators simply introduces additional LECs and we do not investigate this aspect here. We 
assume that k and (3 are the same in full, quenched and partially quenched QCD. Also, and O aa can couple to 
the 77' in QQCD, but the couplings are 1/N C suppressed |15|. and we neglect them. 

The diagrams contributing to /p, , and -Bp (a) are presented in Figs. 0]and[5]respectively. Only diagrams (c) and (d) 
depend on the heavy meson mass shifts. Therefore it is the relative weight between these diagrams and the "tadpole" 
diagrams (a) and (b) which determines the dependence on the heavy meson mass in finite volume effects. 

Although this is the first one-loop calculation for these decay constants and B parameters in finite volume, some 
results in the infinite volume limit already exist in the literature: fp (s) have been calculated at the lowest order in 
full 01 1 Quenched H3. ITEf and partially quenched 0] QCD, and up to first-order corrections in the chiral and 1/Mp 
expansions in full |16| and quenched QCD [l^ ■ The B parameters have been calculated only at lowest order 0, • 
Our results, as presented in Appendix [El agree with all these previous calculations in the appropriate limits. 

B. Phenomenological impact 

We have used the one-loop results in Appendix iBl to investigate the impact of finite volume effects on £. In this 
work, we only intend to estimate the possible size of errors in this quantity, and will leave the actual comparison with 



Notice that the symbol £ as defined in Eq. 1611 is in the traditional notation in B physics, and has nothing to do with the Goldstonc 
field f introduced in Section ITT1 
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— _ ' -X- - . 

i»-0-'* ^MD- - * 

(a) (b) 



{«) (d) 

FIG. 4: Diagrams contributing to the one-loop calculation of decay constants. The thin (thick) solid lines are the heavy-light 
pseudoscalar (vector) mesons. The dashed lines are Goldstone particles, and the crosses are the "double poles" which appear 
in (P)Q^PT. The open squares are the operators defined in Eq. 1621 and the dots are vertices from the HMxPT Lagrangian. 
Diagrams (c) and (d) are for wavefunction renormalisation. 



(a) (b) 



W (d) 

FIG. 5: Diagrams contributing to the one-loop calculation of the B parameters. The open circles are the operators defined in 
Eq. 

lattice data to a future publication. Notice that the one- loop calculation is only valid for <C A x . Nevertheless, we 
present our results for M v up to ~ 800 MeV, where, in principle, higher-order chiral corrections should be included. 
However, finite volume effects are exponentially suppressed at such large pion masses. 

Following the usual procedure in lattice calculations for £, we study two SU(3) breaking ratios 

0- = S -f- (64) 

IB 

and 

Zb = (65) 

in terms of which, 

£ = £/V& (66) 

Furthermore, we define 

K/)fv and (Cb)fv (67) 



to be the contributions from finite volume effects, i.e., those from the volume-dependent part in the one-loop results 
presented in Subsection IIV Al To be more precise, we use the formulae collected in Appendix [B] to calculate the 
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FIG. 6: (£/)fv in full QCD plotted against M v , with L = 1.6 fm. The pion mass M n = 0.35 GeV corresponds to M V L — 2.8, 
and M n — 0.5 GeV corresponds to M^L — 4 in this plot. 



full QCD <£ B ) FV full QCD (£ B ) FV 




FIG. 7: (£b)fv in full QCD plotted against M n , with L = 1.6 fm. The pion mass M v = 0.35 GeV corresponds to M W L — 2.8, 
and M n = 0.5 GeV corresponds to M-^L = 4 in this plot. 



volume corrections with respect to the lowest-order values of fs s {Bb s ) and Jb (Bb), then take the difference between 
the results as an estimate of (£/)fv [(£b)fv]- Since these SU(3) ratios are not very different from unity (at most 
~ 20%), this is a reasonable estimate of these effects. 

Traditionally, many quenched lattice simulations of Bb (s) and were performed using L ~ 1.6 fm. Therefore 
we present our estimate for finite volume effects in QQCD with this box size. For comparison, we adopt the same 
volume for full QCD. As for PQQCD, we work with L = 2.5 fm where most current high-precision simulations are 
carried out ■ Throughout this subsection, we ensure that the condition 

AUL > 2.5 (68) 

holds in all the plots presented here. 

We first discuss the procedure in full and quenched QCD. When studying the light quark mass dependence of 
(£/)fv and (£b)fv, we follow a strategy similar to that in Ref. [331 - That is, we use the Gell-Mann-Okubo formulae 
to express Mk and M v in terms of and M33 (M33, defined in Eq. (|B11|) . is the mass of the fictitious meson 
composed of s and s quarks) 

^ = 2L) (6g) 



and 
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quenched (^)fv quenched (£)fv 




FIG. 8: (C/)fv in QQCD plotted against M n , with L = 1.6 fm and choices of the couplings g and 7. The pion mass = 0.35 
GeV corresponds to M n L = 2.8, and M» =0.5 GeV corresponds to M n L = 4 in this plot. We set a = and M = 700 MeV. 



We investigate the situation where a lattice calculation is performed at the physical strange quark mass (m s ) P hys) but 
the up and down quark mass m is varied. By using (Mx) P hys = 0.498 GeV and (M T ) p h ys = 0.135 GeV we fix, 

(M 33 )p hy s - 2S (m s ) phys = 0.691 GeV, (71) 

as an input parameter in our analysis. Notice that (M33) p h ys is not the mass of a "physical" meson, and the subscript 
just means this mass is estimated by using physical kaon and pion masses. To the same order, we can adopt Eq. <|28[) 
to write 

S s = Ai (M| 3 - Ml) , (72) 

and use (M3 3 ) p hy S , (-Mjr) P hys and physical Mb s — Mb = 0.091 GeV 4] to determine 

Ai = 0.1982 GeV" 1 . (73) 

This determines how 5 S varies with M n . We have also tried to use vanishing pion mass and Md, — Md = 0.104 GeV 
to fix (M 33 ) phys and Ai, and the results presented in this subsection are not sensitive to this variation from the 
values quoted above. 

The results for (£/)fv and (£b)fv for full QCD and QQCD from this analysis are presented in Figs.EHHI with two 
different values for the coupling g (and also 7 in QQCD). Here we stress again that the influence on finite volume 
effects from the presence of A* and S s depends on the size of these couplings, which are not well determined. Inspired 
by the recent CLEO measurement of g in the charm system 0, , and a recent lattice calculation [5(| , we vary g 2 
between 0.3 and 0.5. As for the coupling 7, which is a quenching artifact and has never been determined, we vary 
its value between g and —g. It is clear from these plots that the finite volume effects are generally small in full QCD 
(< 2%), but can be significant in QQCD (~ 3% to ~ 7% for £g) in the range of M^L where lattice simulations are 
normally performed. This is clearly due to the enhanced long-distance effect arising from the "double pole" structure 
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quenched (^b)fv quenched (^b)fv 




quenched (£ b )fv quenched (£ b )fv 




FIG. 9: (£b)fv in QQCD plotted against M n , with L = 1.6 fm and choices of the couplings g and 7. The pion mass M v = 0.35 
GeV corresponds to M W L = 2.8, and M* = 0.5 GeV corresponds to M n L = 4 in this plot. We set a = and M = 700 MeV. 



in (P)QQCD, as first pointed out in Ref. and manifests itself in various places, e.g., nucleon-nucleon potentials 
p?H an d 7T— 7r scattering |rl l38l I3H l52| . 

Although it has been well established that infinite volume chiral corrections are smaller in the B parameters than 
in the decay constants due to the coefficient in front of g 2 in the one-loop results, it is clear from these plots that finite 
volume effects are more salient in £b than in All the quenched lattice calculations for £g have so far concluded that 
this quantity is consistent with unity with typically 3% error. However, we find that the volume effects are already 
at the level of 3—4% when M w = 0.45 GeV in a 1.6 fm box where many quenched simulations were carried out. This 
error depends on both light and heavy quark masses in the simulation, hence is amplified after extrapolating the 
result to the physical quark masses. Also, the fact that volume effects tend towards different directions in full QCD 
and QQCD when M v becomes smaller indicates that quenching errors in these quantities can be larger than those 
estimated in Ref. [l^. Since finite volume effects have not been included in the analysis of lattice calculations of £s 
hitherto, one should be cautious when using the existing quenched results for this quantity in any phenomenological 
work. 

For the analysis in PQQCD, we assume that both the valence and sea strange quark masses are fixed at that of 
the physical strange quark. However, we vary the light sea quark mass to. For this purpose, we define Mss to be the 
mass of the meson composed of two light sea quarks. Therefore, 



M 2 SS 



(Mas)? 



(74) 



/phys V s / m s — physical m s 

Also, we can express the mass shifts S s and 5 sca in terms of meson masses: 

S s = Ax (Af 3 2 3 - Mis) (75) 

and 

<5 SC a = Ai (M j s - Ml) (76) 
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FIG. 10: (£/)fv in PQQCD plotted against Mss (see text for the definition of Mss), with L = 2.5 fm and two different values 
for M w . The pion mass Ml = Mf 3 /4 corresponds to M W L = 4.4 and M^ = M| 3 /8 corresponds to M^L = 3.1. The mass 
Mss = 0.197 GeV corresponds to M S sL = 2.5, and Mss = 0.32 GeV corresponds to M S sL = 4 in this plot. 



by using Eqs. (|2*Tj|) and (fHUI) with the same value of Ai as in Eq. J73J). 

The results for the PQQCD analysis are presented in Figs. EH an d E] The double pole insertions also appear 
in PQQCD and it is clear from these plots that finite volume effects cannot be neglected if one hopes to determine 
£ to the level of a few percent. Especially, in the range of M v and L where current and future lattice simulations 
are performed [47|. they can already be at about 4%, and the dependence on the heavy meson mass is quite strong. 
Therefore they can become comparable to the error presented in the latest review |3lT ] , 

£ = 1.23 ±0.10 (77) 

after quark mass extrapolations. 7 



V. CONCLUSION 



We have investigated finite volume effects in heavy quark systems in the framework of heavy meson chiral perturba- 
tion theory. The primary conclusion is that the scales A* and 5 S , which are heavy- light meson mass splittings arising 
from the breaking of heavy quark spin and light flavour SU (3) symmetries, can significantly reduce the volume effects 
in diagrams involving heavy meson propagators in the loop. The physical picture of this phenomenon is that some 
heavy-light mesons are off-shell in the effective theory, as a consequence of the velocity superselection rule, with the 
virtuality characterised by the mass splittings. The time uncertainty conjugate to this virtuality limits the period 
during which the Goldstone particles can propagate to the boundary. Finite volume effects caused by the propagation 



Finite volume effects presented in this work are, however, correlated with the errors arising from chiral extrapolations. 
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(GeV) 




FIG. 11: (£b)fv in PQQCD plotted against Mss (see text for the definition of Mss), with L — 2.5 fm and two different values 
for M n . The pion mass Ml = M| 3 /4 corresponds to M W L = 4.4 and Ml = M| 3 /8 corresponds to M W L = 3.1. The mass 
Mss = 0.197 GeV corresponds to M S sL = 2.5, and M S s = 0.32 GeV corresponds to M S sL = 4 in this plot. 



of the Goldstone particles naturally affect the light quark mass extrapolation/interpolation in a lattice calculation. 
On top of this, our work implies that they also influence the heavy quark mass extrapolation/interpolation, since 
the scale A* varies significantly with the heavy meson mass. The strength of this influence is process-dependent, 
determined also by the relative weight between diagrams with and without heavy meson propagators in the loop. 
The volume effects can be amplified by both heavy and light quark mass extrapolations. Therefore it is important to 
perform calculations to identify these effects in phenomenologically interesting quantities. 

We have presented an explicit calculation in finite volume HMxPT for the B parameters in neutral B meson mixing 
and heavy-light decay constants, in full, quenched, and partially quenched QCD. We have used these results to estimate 
the impact of finite volume effects in the SU(3) ratio £, which is an important input in determining the magnitude 
of the CKM matrix element Vtd- Within the parameter space where most quenched lattice calculations have been 
performed, we find that, although this impact is quite small (<~ 2%) in full QCD, it can be significant in QQCD. This 
is due to the enhanced long-distance effects arising from the double pole structure. This error will be amplified by the 
quark mass extrapolations and hence can exceed the currently quoted systematic effects. Furthermore, finite volume 
effects tend towards different directions in full and quenched QCD for decreasing M n . This means that quenching 
errors in £ may be significantly larger than what was estimated before. Therefore one has to be cautious when using 
the existing quenched lattice QCD results for £ in phcnomenological work. In PQQCD, our results indicate that finite 
volume effects are typically between 3% and 5% in the data range of future high-precision simulations, and they can 
be significantly amplified in the procedure of quark mass extrapolations. This means that they are not negligible in 
future lattice calculations of £. 
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APPENDIX A: INTEGRALS AND SUMS 

We have regularised ultra-violet divergences that appear in loop integrals using dimensional regularisation, and 
subtracted the term 



A 



4-d 



7£ + log(47r) + 1. 



(Al) 



The integrals appearing in the full QCD calculation are defined by 

d d k 1 



I- X (m) = fi 



16tH 



(2it) d k 2 — m? + ie 



A — log 



(A2) 



where 

F x (m,A) 

with 



16tt 2 



d 



d d k 



OA J (27r) d (k 2 ~ m 2 + ie)(v ■ k - A + ie) 
d 

3— F x (m,A), 



2A 2 



/10A 2 4m 2 A 2(A 2 -m 2 ) /A 
m 2 I A + ( — ; — ) A + — -mR — 



\ 9 



(A3) 



(A4) 



R(x) = \/ x 2 — 1 log 



y/x 2 — 1 + ie 



(A5) 



x + \/ x 2 — 1 + ie I 

and n is the renormalisation scale. For the quenched and partially quenched calculations, we also need the integrals 



d d k 



1 



= dl x {m) 

(2ir) d (k 2 — m? + ie) 2 dm 2 ' 



(A6) 



and 



HUm,A) = {gP v - v p v v )^~ d 



d d k 



kpky 



d_ 

dA J J2nj d (k 2 - m 2 + ie) 2 (v ■ k - A + ie) 
^H x (m,A). 



dm 2 



(A7) 



In a cubic spatial box of size L in four dimensions with periodic boundary condition, so that the three-momentum is 
quantised as in Eq. I)4ip. one obtains the sums (after subtracting the ultra-violet divergences) 



m z + ie 



I(m) + Jfv("i), 



(A8) 



18 



and 



H(m,A) = (gP" - v p 



nf 1 W d f dk ° 



hph u 



2tt (k 2 — to 2 + ie) (v ■ k — A + ie) 



ff (to, A) + H F v(m, A) (A9) 



for the full QCD calculation, where 



and 



I(m) = h(m)\- x=0 , 



H(m)=H x (m,A)\ x=0 , 
are the infinite volume limits of X and Ti, and (n = \n\) 



hv{m) 



47r 2 vz — ' nV 

n#0 



4^ 2^ Y 2n~L 



run ( 1 \ _„„,, , 3 
e nmL x < 1 + 
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8nmL 128(titoX) 2 



O 



(A10) 



(All) 



n 3 N 



nmL 



(A12) 



is the finite volume correction to /(to). The function Hpy is the finite volume correction to H(m, A) and can be 
obtained via 



Hp\?(m, A) = i 



(to 2 - A 2 )K FV (m, A) - 2AJ fv (to, A) + H FV {m) 



where Jpy(m, A) and K-pv{m, A) are defined in Eqs. (|44|l and (|55|) . 
For QQCD and PQQCD calculations, one also needs 



L 3 2n 



1 



dl(m) dlpy(m) 



(fc 2 — to 2 + ze) 2 dm 2 



dm 2 



(A13) 



(A14) 



and 



W (to, A) . - 



kph v 



L 3 J ^ J 27r (k 2 - m 2 + ze) 2 (t> • - A + ie) 

k 



dH(m,A) dH F v(m,A) 



dm 2 



dm 2 



(A15) 



APPENDIX B: ONE-LOOP RESULTS FOR DECAY CONSTANTS AND B PARAMETERS 



In this appendix, we collect the results for one-loop corrections to fp (s) y/Mp (a) and Bp (s) . For convenience, we 
introduce 



and 



C±{M GP ,x) =l(M G p)±g 2 H(M GP ,x), 



0% (M GP , x) = T>' (M gp ) ± g 2 W (M GP ,x), 



(Bl) 



(B2) 



where the functions I(m), H(m,x), I v (to) and Tt v (to, a;) are defined in Eqs. (|A8|) . 1A9|I . 1|A14(I and (|A15|) . respec- 
tively. 
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In full QCD, we find 

/pVMp = nl 1 



12/- 

i 



3/ 2 



Bp = 



Br - 



3£ 
2k 2 

3/3 
2^ 



6.f 2 
2i 

3T 2 " 



9 C_(M ff ,A,) + 6 C-(M K ,A* + 6 S )+C-{M V ,A*) 
3 C-(M K , A* — <5 S ) + C-(M V , A. 
3C + (M 7r ,A*)+C+(M r? ,A : 



C+tM^A*) 



In QQCD, wc find 



}■ 



|c_(M 33 , A«) + f aM ^ M ° ) ^ (M :V ,. A ) r2 // ,H(,.U,,.Aj 



Bf 



2 K 2 \ /2 
2 K 2\ /2 



(l-|)c+(M W) A,)- 
(l-|)c+(M33,A.) 



' ' 1 ' " M ° 1 (M W) A,) + 2 37 W(M T , A,; 



where 



M 33 



2M%-M*. 



In PQQCD, we find 
fpy/Afr = J 1 



W 2 



12 C_(M y5 , A, + 4oa) + 6 C_(M r , A* + <5 sea + <5 S ) 



M| 3 - Ml 25 2 VSs 



M 2 -M 2 3 - 2^ + 2^/ 2 



(B3) 
(B4) 

(B5) 
(B6) 

(B7) 
(B8) 
(B9) 
(BIO) 

(Bll) 



(B12) 



^ s (M| 3 -M 2 -24 Ss ) v . 
+ 6 M2-M| 3 + 4 s + 24 Ss C - (M - A * ) 



fp s VM~ F 



n{ 1 



6/ 2 



3 C-(M VSs , A, + <5 soa + <5 S - (5 S ) + 6 C_(M Cs , A, + (5, 
Mf 3 - M 2 + 25= s 



-C_(M 33 ,A*) + 2 



M| 3 -M 2 + 2^ + 4^ 

^ Ss (M| 3 -M 2 + 2<% g ) 



M 2 - Mj 3 - 251 s + Kss 
Ml - M| 3 - 2<S* s - 452 Sa 



C_(M X ,A*) 



+ 6 



Ml - M| 3 - - 4^ Ss 



(B13) 



B F 



i^/i__L 

2k 2 1 6/ 2 



6C+(M x ,A*)-3 



M 2 



M| 3 + 2*2 5s 



M 2 - M| 3 + SI 



VS "T ^^ygg 



C+(M W ,A.) 



+ 



M 2 



M 2 - M| 3 + # s 



2 ^s + ^vss x C+ (M X ,A*) 



28 2 



VSs 



Sl s (M? 



+ 6 M| - Mf 3 + £ 2 /fl + 2# s > (M - A * } 



M 



(B14) 
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C+(M 33 , A» 



(B15) 



where 



M 2 T = B (m + m s ) = M 2 S 




M vSs = B (m s + m s ), 

M 2 Cs = B (m s +rh)= M* s + M% - M$, 

S VSs = M 33 ~ M VSs, 



and Ml 



i {Ml + 2M 3 2 3 - 25 2 vs 4S 2 VSs ) . 



(B16) 



It is straightforward to show that the PQQCD results reproduce those for full QCD in the limit to = to and to s = to, 
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